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Abstract 

Using the language of non-relativistic effective Lagrangians, we formulate a systematic framework 
for the calculation of resonance matrix elements in lattice QCD. The generalization of the Liischer- 
Lellouch formula for these matrix elements is derived. We further discuss in detail the procedure of the 
analytic continuation of the resonance matrix elements into the complex energy plane and investigate 
the infinite- volume limit. 
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1 Introduction 



The calculation of matrix elements involving unstable states has already been addressed in 
lattice QCD. As examples, we mention the recent papers [1-3], which deal with the electromag- 
netic form factor of the p-meson, as well as the electromagnetic and axial-vector form factors 
of the A-resonance and the NA'y transition vertex. Electromagnetic and axial transition form 
factors for the Roper resonance have also been studied [4]. Moreover, we expect that the num- 
ber of such investigations will substantially grow in the nearest future due to a growing interest 
in the study of the excited states. 

Even if one argues that the quark (pion) masses in the above lattice simulations are large, 
so that all resonances are in fact stable particles, various conceptual questions arise: 

i) It is clear that we are ultimately interested in simulations carried out at the physical quark 
masses. Is it possible (at least in principle) to tune the quark mass continuously until it 
reaches the physical value? 

ii) In the continuum field theory, any matrix element with resonance states is defined through 
an analytic continuation of the three-point Green function into the complex plane P 2 — > 
sr, where P denotes the pertinent four-momentum and sr is the resonance pole position 
in the complex plane (its real and imaginary parts are related to the mass and the width 
of a resonance). What is the analog of this procedure in lattice field theory? 

iii) Once this procedure is defined, what is the volume dependence of the measured form 
factors? 

In this paper, we address these questions in detail. In order to formulate the problem in 
a more transparent manner, let us first define what is meant by resonance matrix elements in 
the continuum field theory and on the lattice. We start with the continuum field theory and, 
for simplicity, concentrate on the scalar case. Consider an arbitrary (local or non-local) scalar 
operator 0(x) which has the internal quantum numbers of a given resonance. The statement 
that a resonance is present is equivalent to the claim that the two-point function 

D(P 2 ) = % J d 4 xe lPx {0\TO(x)O\0)\0) (1) 

has a pole in the complex variable P 2 on the lower half of the second Riemann sheet at P 2 = sr: 

lim D(P 2 ) = + regular terms at P 2 —> sr. (2) 

P 2 ^s R Sr - P 2 

The real and imaginary parts of sr are related to the resonance mass tur and the width Tr, 
according to ResR = m 2 R — T 2 R /A, Ims R = — itirTr. 

In order to define resonance matrix elements 1 , say, of the electromagnetic current J M , we 
consider the following three-point function: 

^(P,Q) = i 2 [ d A xd*ye iPx - iQy (0\TO(x)J^0)O\y)\0) . (3) 



1 The following discussion is a straightforward adaptation of the procedure which has been used to define the 
matrix elements in the case of stable composite objects, sec, e.g., Rcfs. [5,6]. 
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The form factor of a resonance is then defined as 

i(P + Q) ll F((P-Q) 2 ) = (P\J,(0)\Q} = lim B R 1/2 (s R - P 2 )^(P, Q)(s R - Q 2 )B R 1/2 , (4) 

P 2 ,Q 2 ^s R 

where Br is the residue at the resonance pole, see Eq. (2). Note that the matrix element 
displayed in Eq. (4) should be understood as a mere notation: in the spectrum, there exists no 
isolated resonance state with a definite momentum. Moreover, as it is clear from Eq. (4), this 
definition of the resonance matrix elements necessarily implies an analytic continuation into the 
complex plane. We would like to stress that we are not aware of any consistent field-theoretical 
prescription, where the analytic continuation would not be employed. 

Let us now ask the question, how such resonance matrix elements could be evaluated on 
the lattice (at least, in principle). As it is well known, a resonance does not appear as an 
isolated energy level. There exist alternative approaches to the problem of extracting resonance 
characteristics (the mass and the width) from the measured quantities on the Euclidean lattice. 
In this paper, we work within Liischer's finite- volume framework [7] 2 . In order to calculate the 
matrix element on the lattice, one usually considers the following three-point function 

V,(P,t';Q,t) = (0|TO P (O^(0)Oj,(t)|0>, (5) 

where 

Qp(0 = ^ Px O(x, , 0{,(i) = e iQy O f (y, t) . (6) 
x y 

In addition, we define: 

D(P,t) = (0|rOp(t)Oj,(0)|0>. (7) 

The matrix element of the electromagnetic current between the ground-state vectors in a channel 
with the quantum numbers of the operator 0(x), moving with the 3- momenta P and Q, 
respectively, is given by 

W )| W | W ) = f Mit(P,/;ftfljZ^^IZ. (8 ) 

t— >— OO 1 

2 At present, Liischer's approach [7] has been widely used to obtain scattering phase shifts from the energy 
spectrum in a finite volume. The resonance position can be then established by using the measured phase shift. 
The procedure can be directly generalized to the case of multi-channel scattering [8-10]. Moreover, in Ref. [11] 
it has been argued that the use of the physical input based on unitarized Chiral Perturbation Theory may 
facilitate the extraction of the resonance poles from the lattice data (the method has been subsequently applied 
to different physical problems in Refs. [12, 13]). Recently, a generalization of Liischer's approach in the presence 
of 3-particle intermediate states has been proposed [14]. Other approaches to the determination of resonance 
pole positions imply the study of the two-point function at finite times [15,16], as well as reconstructing the 
spectral density by using the maximal entropy method [17]. The application of different approaches to the 
extraction of the resonance properties from the lattice data has been carried out recently in Ref. [18]. Last but 
not least, the finite volume approach has been applied to study the two-particle decay matrix elements on the 
lattice [19-21], including the case of multiple channels [22]. 
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Using the generalized eigenvalue equation method, the matrix elements between the excited 
state vectors \E n (P)) can be also defined in a similar manner. 

If the ground state of a system corresponds to a stable particle, then Eq. (8) indeed yields 
the form factor of a stable particle in the infinite- volume limit, which in this case is well defined. 
However, the situation in case of resonances is conceptually different. The easiest way to see 
this is to note that in the infinite- volume limit the energy of any state \E n (P)) tends to the two- 
particle threshold energy. In other words, any given energy level decays into the free particle 
levels in the limit L — > oo (here, L denotes the size of a spatial box). Moreover, as shown 
in Ref. [23] (in case of 1 + 1 dimensions), the matrix elements measured for any given level 
follow a similar pattern. For example, the magnetic moment tends to the sum of the magnetic 
moments of the free particles in the limit L — > oo. Obviously, this is not the result that we 
wish to extract from lattice data. 

As mentioned above, using Luscher's approach, it is possible to determine the resonance 
pole position in the complex plane from the measured finite-volume (real) energy spectrum. 
This position stays put (up to exponentially suppressed corrections in L) in the limit L — )■ oo, 
despite the fact that all individual levels collapse towards threshold in this limit. The aim of 
the present paper is to formulate a similar approach for the matrix elements, and to ensure 
that the matrix elements that are extracted with the help of such a procedure coincide with 
the infinite- volume matrix elements, e.g., given in Eq. (4), up to exponentially suppressed 
corrections. 

The goal, stated above, will be achieved by a systematic use of non-relativistic effective 
field theory (EFT) in a finite volume. In particular, we shall calculate the quantity in Eq. (8), 
which can be measured on the lattice, within non-relativistic EFT, and shall identify a piece 
in this expression, whose infinite-volume limit coincides with the resonance matrix element in 
the infinite volume we are looking for. 

The paper is organized as follows: In section 2 we formulate a covariant non-relativistic EFT 
in a moving frame and re-derive the Gottlieb-Rummukainen [24] formula within this approach. 
The extraction of a resonance pole position is discussed in detail. In section 3 we give a short 
re-derivation of the Luscher-Lellouch formula [19], as another application of the non-relativistic 
EFT methods. Further, in section 4 we evaluate the vertex function in the non-relativistic EFT. 
The infinite- volume limit of different terms in the expression of the vertex function is analyzed 
in detail in section 5, where particular attention is paid to the so-called fixed singularities that 
emerge in a result of analytic continuation of Luscher's zeta-function into the complex plane 
in 3 + 1 dimensions. The prescription for calculating the resonance matrix elements is given in 
section 6. Section 7 contains our conclusions. 

2 Extraction of the resonance poles in moving frames 

The initial and final states in a form factor have non-zero momenta. For this reason, one has to 
formulate a procedure for extracting resonance pole positions in moving frames. Within poten- 
tial quantum mechanics, this has been done in Refs. [24], see also Ref. [25,26] for the general- 
ization to the non-equal mass case. Refs. [20, 27] address the same problem in a field-theoretical 
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setting. Finally, in Ref . [28] , a full group-theoretical analysis of the resulting equation has been 
performed, including the case of particles with spin. Below, we shall briefly re-derive this result 
within the non-relativistic EFT along the lines similar to Refs. [29, 30], where the treatment was 
restricted to the rest frame. At a later stage, the same approach will be used for the calculation 
of the matrix elements. 

In the treatment of the moving frames it is very convenient to use the covariant form of 
the non-relativistic EFT which has been introduced in Ref. [31] and was discussed in detail 
in Ref. [32]. Assume, for simplicity, that we deal with two elementary scalar fields $12 with 
masses m^, respectively. The Lagrangian is given in the following form: 

i=l,2 

+ d (($ t 1 )^($ 2 ) At $ 1 $ 2 - m 1 m 2 ${$ 1 $l$ 2 + h.c.) 

+ C 2 (^ 2 r ~ ($l) M $5)(($i)„$2 - *i($2)„) + • • ■ , (9) 

where $j, i — 1,2 denote the non-relativistic field operators, Wi = \/mf + A are the energies 
of the particles (here, A = V 2 ), and 

(**) M = (7> 4 ) M *i , = {V\)M , (P*) M = (Wi, -tV) , = (Wi,iV). (10) 

Further, the ellipses stand for terms containing at least four space derivatives. To set up 
the power-counting rules we introduce, as in Refs. [31,32], a generic small parameter e and 
count each 3- momentum as = 0(e), whereas the masses are counted as = 0(1). The 
Lagrangian given in Eq. (9) contains all allowed explicitly Lorentz-invariant terms 3 up-to-and- 
including 0(e 2 ), and the omitted terms are of order e 4 . 

The non-relativistic couplings Cq, C±, O2, ■ • • , which are present in the Lagrangian, are di- 
rectly related to the effective-range expansion parameters for 1 + 2 — > 1 + 2 elastic scattering 
(scattering length, effective range, etc), see Refs. [31,32]. We would like to remind the reader 
here that the theory described by the Lagrangian given in Eq. (9) conserves particle number, 
so it can be applied in the elastic region only. 

The Feynman rules, which are produced by the Lagrangian (9), should be amended by 
a prescription which states that the integrand in each Feynman integral is expanded in 3- 
momenta, each term is integrated by using dimensional regularization and the result is summed 

3 Note that in the conventional non-relativistic theory the number of the allowed terms at a given order in 
e is much larger, because these terms are not restricted by the requirement of Lorentz-invariance. At the end, 
however, matching to the relativistic amplitude should be performed that effectively imposes such constraints on 
the low-energy couplings, because the number of physically independent low-energy parameters in the relativistic 
amplitude is smaller. In this way, the constraints are imposed in a perturbative manner, order by order in e. 
On the contrary, in our approach, we impose the requirement of the Lorentz invariance from the beginning and 
avoid the introduction of the constraints at all. The key property which allows us to do this is that in our 
approach (unlike the conventional framework) non-relativistic loops are Lorentz-invariant by itself, so it suffices 
to impose Lorentz-invariance at tree level only. For more details, we refer the reader to Ref. [32]. The method 
of matching to the relativistic theory was already used in the construction of the heavy-baryon chiral effective 
Lagrangian in Ref. [33] . 
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up again [31,32]. Below, we shall consider the theory in a finite volume. It is easy to see that, 
for consistency, one should apply the same prescription, replacing the dimensionally regularized 
integrals by sums over discrete momenta. In particular, one has to discard everywhere discrete 
sums over polynomials in momenta, in accordance with the similar infinite-volume prescription 
in the dimensionally regularized theory. 

Let us start in the infinite volume. Using the above Feynman rules, it is straightforward to 
ensure that the scattering T-matrix in the infinite volume in an arbitrary moving frame obeys 
the Lippmann-Schwinger (LS) equation: 

T{pi, p 2 ; qi, q 2 ) = -V (pi, p 2 ; qi, q 2 ) 



(27r) d 2wi(ki) (27r) d 2u; 2 (k 2 ) 



x (2tt) 6 (pi + p 2 - ki - k 2 ) — — — — — — -, (11, 

wi(ki) + w 2 [k 2 ) ~ wi(pi) - w 2 {p 2 ) - i0 



where Wi(Y) = \J mf + l 2 and the potential is given by the matrix element of the interaction 
Hamiltonian, which is derived from the Lagrangian (9) by the canonical procedure, between 
the two-particle states 

(pi,p 2 |#i|qi,q 2 ) = (2tt) 3 (5 3 (pi + p 2 -qi -q 2 )^(pi,p 2 ;qi,q2)- (12) 

Note that we have used dimensional regularization in Eq. (11). The parameter d denotes the 
number of space dimensions (at the end of calculations, d — > 3). 

By construction, the potential V is a Lorentz-invariant low-energy polynomial that depends 
only on scalar products of the 4-momenta. The first few terms in the expansion are given by 

-V^Pi, P2; qi, q 2 ) = C + Ci(pip 2 + <?i<? 2 - 1m x m 2 ) + C 2 (p 2 - Pi)(<? 2 - <?i) + 0(e 4 ) , (13) 

where, e.g., pip 2 = Wi(pi)w 2 (p 2 ) — Pip 2 , etc. In general, defining the center-of-mass (CM) and 
relative momenta, according to 

P = Pi+P2, P = V2P1 ~ , Hi,2 = o 1 ± ™ ' P ~ 



2 V P 2 J 4:P 2 
Q = qi + q2, q = V 2 qi ~ \h.q% , Mi,2 = „ 1 ± 7v? ' 9 = 77r? A 14 ) 



2 V Q 2 J 4Q 2 

where X(x,y,z) denotes the Kallen triangle function, it can be seen that V is a low-energy 
polynomial of six independent Lorentz-invariant arguments p 2 , q 2 , pq, PQ, Pq, pQ. The original 
arguments p±p 2 , Piq\, Pi 92) P2qi, p 2 9 2 , qiQ2 can be expressed through linear combinations of these 
arguments with coefficients, which themselves are low-energy polynomials. 

Consider now the partial-wave expansion of the potential. To this end, we define the mo- 
menta boosted to the CM frame (note that the boost velocity is different in the initial and the 
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final states, because the potential is generally off the energy shell): 

p = p + PI (7 — 1)— — 7t!jpj I , p = 7 p - 7 w— 

p; = (Vp*o), « = E 7 = (i-^T 1/2 , 



q =q + Q^(7-l)^-7«|Q|J 



/ / /qQ 

9o = 7 ?o - 7 u TqT 



q; = (^,o), «' = M, y = (i - K)T 1/2 ■ (is) 



Qo 

Taking into account the fact that P = Q in the "lab frame," it is straightforward to show that 
pq = p* q * + O((P -Q ) 2 ) = -p*cC + O((P°-Q ) 2 ), 



PQ = =P*Q* + 0((P -Q ) 2 ) = VpZy^ + 0((Po-Qo) 2 ), 
Pq = p* q * + o(P Q -Q ) = + 0(P -Qo), 

P Q = p*Q* + O(P -Q )=0 + O(P -Q ). (16) 

In addition, p 2 and q 2 can be expressed in terms of P 2 and Q 2 , respectively. This means that, 
up to terms that vanish as Po — ^ Qo-, the potential can be rewritten in the following form 

-^(Pi,p 2 ;qi,q 2 ) = -47r5>(|p*|, |q |)^ m ( P *)X m (q ) + O(P - Q ) . (17) 

Im 

Here, the function v\ can be chosen to be real and symmetric with respect to its arguments, 
i.e., Eq. (17) describes a Hermitean potential. The quantity 3^ m (p) is defined as 3^m(p) — 
|p|'yj m (p), where Yi m are the usual spherical harmonics. The terms that vanish as Po — > Qq 
can be omitted from now on. The justification for this is the fact that the parameters in the 
potential are determined by matching to the physical S'-matrix elements (on shell), order by 
order in the low-energy expansion. The omitted terms do not contribute either at tree level or in 
loops (the latter because the regular momentum integrals vanish in dimensional regularization). 
Consequently, one may consistently set these terms equal to zero from the beginning. 
Performing now the partial-wave expansion in the amplitude 

T( Pl ,p 2 ; qi ,q 2 ) =47r^tKIP*l > |q*l)^m(p*W m (q*), (18) 

im 

substituting this expansion into the LS equation (11), and using the properties of dimensional 
regularization, on the energy shell |p*| = |q*| = A 1 ^ 2 (s, mf, m 2 ) / (2y/s) we get 

U{s) = -in(s) - ^(s)|p*| 2 'G( s )tKs) , (19) 
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where the obvious shorthand notations for the on-shell quantities vi(s) = i>{(|p*|, |p*|) and 
ti(s) = ti(\p*\, |p*|) are used. The quantity G(s) is given by [31,32]: 

f gki d d k 2 (2ir) d 5 d (P - ki - k 2 ) i\p*\ 

[S) J (2 7 r) rf 2 U ; 1 (k 1 )(27r) rf 2 M ; 2 (k 2 ) w^) + w 2 (k 2 ) - P - iO 8n^' 1 } 

Further, unitarity gives: 

where Si(s) is the scattering phase. 

The transition to the finite volume is performed in the "lab frame". The momenta are 
discretized according to 

^ = j- n, , n, G Z 3 . (22) 

The partial-wave expansion of the potential does not change. However, since the introduction of 
a cubic box breaks rotational symmetry, the partial- wave expansion of the scattering amplitude 
has to be modified: 

T( P i,p 2 ;qi,q 2 ) = (4tt) £ t lm>Vm ,(\p*\, |q*|; P)^ m (p*)^ m ,(q*) . (23) 

lm,l'm' 

Substituting this expression into the Lippmann-Schwinger equation, on the energy shell we 
obtain: 

tl m ,l'm'(s; P) = -5lmJ'm'Vl{s) - ill ^ V\ {s)X [m ^ m „ (s, P)Vm",«'m' ( s ! P ) > ( 24 ) 

i"m" 

where 

x f s p)-iv y;jk*)y,, m ,(k*) , . 

Wl ' j " ^ 3 2 ^i( k i) 2 ^( p - ki)K(kO + W2 (P - kO - P ) • 1 j 

Next, we use the identity [32] 

1 1 1 



4u>iu> 2 (u>i + w 2 - P ) 2P (kP) 2 

K P 2 IP J 

+ -J-f ! i ! V(26) 

4w!tt; 2 \tu! + u? 2 + P wi-w 2 + Pq w 2 - w 1 + P J 

where k = ki + /iiP. One can straightforwardly check that the term in the brackets does 
not become singular in the physical region. Using the regular summation theorem [34], one 
may then replace the sum over ki in this term by the integral. Further, to be consistent with 
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our prescription for the calculation of the Feynman integrals in dimensional regularization, 
one should put these integrals to zero. After this, the expression for Xim t i> m '(s, P) takes the 
following form: 

Xlm,l>m>(s : P) - — ^ (27) 

k=k 1 +A t iP k 2 - i — f- - (p*) 2 



In order to transform this equation further, let us define the parallel and perpendicular com- 
ponents of the three vectors with respect to the CM momentum P. In particular, one may 
write k* = (Jfe[,kl), where fcjj = (7*) _1 Jfc||, = k ± and 7* = (1 - (v*) 2 )~ 1/2 , v* = \P\/E* = 
|P|/(^i(k*) + 1^2 (k*)). Consequently, on the energy shell E* = Pq we obtain: k* = r = 
(7~ 1 fc|j, kj_) with 7 = (1 — P 2 /P 2 ) _1 / 2 . Up to exponentially suppressed terms, Eq. (27) now 
takes the form 

(p*y+ l '+ 1 _, 

M W 3 



■J'm' 

' j=\i-i'\ s =-i 



C lm , js , M = {-r't*-i+*y/(2l + l)(2j + 1)(2J' + 1) ( l m I J m )( l J Q [), {2i 



where 



and 



-|P, (29) 



Zt L (l; s) = , P, = {r = M 3 |r | |=7" 1 H-/ii|d|), r ± = n ± , n G Z 3 } . (30) 

Note that Zf m (l; s) is a function of s and not merely rj 2 , as in the rest frame. This happens 
because the kinematical factor 7 depends on s. 

The finite-volume spectrum is determined by the pole positions of the scattering matrix. 
The poles emerge when the determinant of the system of linear equations (24) vanishes. Taking 
into account Eqs. (21,28), the equation determining the the energy spectrum can be written in 
the following form: 

det(s w 5 mm/ - tan 5i(s)Mim,i'm'(s; P)j =0. (31) 

This is Liischer's equation in a moving frame, or the Gottlieb-Rummukainen formula (see 
Refs. [20,24,25,27]). It can be also shown that, in the large- L limit, the equations obtained 
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in Ref. [35] reduces to Eq. (31), if all partial waves, except the S-wave, are neglected. Using 
discrete symmetries, the system of linear equations (31), that couples all partial waves, can be 
partially diagonalized. We do not, however, address this problem here. A full-fledged group- 
theoretical analysis of the Gottlieb- Rummukainen formula with the inclusion of the spin of the 
particles forms the subject of a separate investigation [28]. 

The equation (31) enables one to extract the scattering phase shift from the measured energy 
spectrum on the lattice. In order to extract a resonance pole position in the complex plane 
from the phase, additional effort is needed. For example, one could assume that the effective 
range expansion is valid up to the resonance energy. This assumption works well, e.g., for the 
physical A-resonance. The effective-range expansion for the scattering phase shift is written 
as: 

P 2l+1 cot *,(«) = -- + \ r lP 2 + 0(p 4 ) , p 2 = A( ^ X) . (32) 
a\ 2 4s 

This means that the lattice data allow one to determine the scattering length a/, the effective 
range rj, etc. The pole position p R (on the second sheet) is then determined by solving an 
algebraic equation with known coefficients: 

pl +l cot 6i(s R ) = -~ + \ npl + ■■■ = -ip* +1 . (33) 
ai I 

It should be stressed that, in order to justify the application of this procedure, the data should 
cover the energy range where the resonance mass is located. There exist alternative strategies, 
which may be applied, if the use of the effective-range expansion is questionable. However, 
the present paper is mainly focused on the study of the resonance matrix elements. In order 
to make the conceptual discussion of this issue as transparent as possible, below we restrict 
ourselves to the situation where the effective-range expansion can be used without problems. 

3 Liischer-Lellouch formula for the scalar form factor 
from the non-relativistic EFT 

Before investigating the resonance matrix elements, we consider the simpler problem for matrix 
elements of stable states and re-derive the Liischer-Lellouch formula [19] in an arbitrary moving 
frame within the non-relativistic EFT. To ease notations, we treat the equal mass case m\ = 
m 2 = m here, albeit the formalism can be straightforwardly generalized to the unequal-mass 
case 4 . As an example, we consider the (scalar) form factor in the time-like region. In order 
to study the form factor, the non-relativistic Lagrangian in Eq. (9) should be equipped by the 
part that describes the interaction with the external field A(x). This part of the Lagrangians 
takes the form 



£ A = eA(x)j(x) = eA{ <f>\<f>\ + D x 



$ t 1 ($^ + 2 ($t r( $t )/j + (<i) t ) M $ t 



+ • • \ + h.c. , (34) 



4 We have in mind, e.g., the calculation of the pion form factor. 
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+ + + + 



Figure 1: Diagrams contributing to the matrix element on the l.h.s of Eq. (37) in perturbation 
theory. The first diagram corresponds to the free propagation of the fields $1,2 ■ 



where (cf with Eq. (10)) 

(^::. = (vlY---(vX---®l (35) 

and the low-energy constants e,D\,--- describe the coupling of the field A(x) to $12 (note 
that a similar approach to the electroweak matrix elements in the two-nucleon sector of QCD 
was adopted in Ref. [36]). 



Define now the operators 

fL/2 

O(x ;P,k) = / d 3 xd 3 ye-^ P(x+y) - ?k(x " y) $i(xo,x)$ 2 (xo,y) 

J-L/2 
fL/2 

O*(x ;P,k) = / d 3 xd 3 ye^ P(x+y)+ik(x_y) $l(xo,x)$5(xo,y), 

J-L/2 



-L/2 

^ 27T 27T / 1 \ „a 

P = — m, k=— n+-m , m, n G Z 3 , (36) 

Lj ±j \ Z J 

and consider the following matrix element in Euclidean space for xo > yo- 

(0\O(x ; P, k)O\y ; P, k)|0) = £ |(0|O(0; P, k)|£ n (P))| 2 e - E ^°-^ , (37) 

n 

where the E n = E n (P) denote the energy eigenvalues for the eigenstates with total momentum 
P. 

Note that in the non-relativistic EFT the above matrix element can be calculated in per- 
turbation theory. The pertinent diagrams are shown in Fig. 1. Using the Euclidean-space 
propagator in the non-relativistic EFT 

r rjn^ 1 P *PoOo-yo)+«p(x-y) , 

WHx)*t(ym = j^ T3 Z 2a{p){a{p) + ipa y = (38) 

for this matrix element we get: 

(0\O(x ; P, k)O\y 0] P, k) |0) = L 3 j ^ e l 



p 



3 / dP JP (x -y ) 



-I 



L 3 



4 Wl (k) W2 (k)(p - i( Wl (k) + w 2 (k))) (4 Wl (k)w 2 (k)y(p - i( Wl (k) + w 2 (k))y 
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(39) 



where 

Wi(k) = ^m2+(! + k) , W2 (k) = ^m 2 +(|-k) , (40) 

and T is the forward scattering amplitude of the particles 1 and 2 in the moving frame (see 
Fig. 1): 



T C » + C o L3 E 4t(;i(1W1)(wi(1) + aj2(1)+l p o) + 



C 'o + C , o 87r Jy i7 ^o d o(l^) + -- 



-1/2 



-OPo+P 2 ), 7=U + 7*) , P* = ^-m 2 , V = ^, (41) 



P 2 / ' J V 4 ' ' 2tt 

where we have used Eqs. (25,28), and where we have retained only the S-wave contribution 
in the scattering matrix in order to simplify the discussion of the scalar form factor. Using 
Eqs. (13,21), the tree-level and bubble diagrams in Fig. 1 can be summed up to all orders. The 
result on the energy shell is given by 

tan0 d (s) = - . '\ , (42) 



p* cot 5(s) +p* cot 4> d (s) ' Z^ (l;s) 

where S(s) = S (s) denotes the S-wave phase shift. 

The eigenvalues are determined from the Gottlieb- Rummukainen equation (see section 2): 

S( s ) = -0 d (s) + irn , s = s n , P fixed. (43) 

The quantity T defined by Eq. (42) has poles at real values of s = s n , i.e., at P = Pon = iE n 
where E n = E n (P) = s n + P 2 . In the vicinity of this pole, the quantity T behaves as: 



327rsin S(s n ) 

T — > — — — — — — h regular terms , (44) 

6>(s) + (0 d (s„))' E n P n + zP 1 ; 

where the derivative is taken with respect to the variable p*. Substituting now this expression 
into Eq. (39), performing the integral over P and taking into account the fact that the "free" 
poles at Pq = i(w\(k) + u>2(k)) cancel in the integrand, the final expression for the matrix 
element in Eq. (39) for xq — yo > reads: 

,2 



(0|O( l0 ;P,k)Ot (!/0 ;P, k )|0) = L'Y. ^TJul ,, -„ 

3— yo) 

(45) 



5'(s) + ((f> d (s n )y E n 

e -E n (x -y ) 



(A Wl (k)w 2 (k)y(E n - Wl (k) - u, 2 (k))2 
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Figure 2: Diagrams contributing to the vertex function. The wiggly line corresponds to the 
external field A(x). 

Comparing this expression with Eq. (37), one reads off: 

1/2 



|<0| O (0 ; P, k )| £ „(P»| = ^( ^^yi f 

1 1 

X A Wl {k)w 2 {k) \E n - Wl (k) - w 2 (k)\ ■ ( 6) 

Next, we turn to the determination of the form factor in the time- like region. To this end, we 
have to consider the amplitude of pair creation from the vacuum in the presence of an external 
field A(x), at the first order in the coupling e. This matrix element is described by 

(0\O(x , P, k) £ A (0)|0) = eA(0)F[k, P; x ], x > . (47) 

We evaluate the quantity F in perturbation theory. The pertinent diagrams are shown in Fig. 2. 
Summing up all bubbles yields: 

Mk p. _ i = p (f) f dP, e^> f cot 5{s) 

1 ' ' 0J y, J 2m4w 1 (k)w 2 {k)(P -t(w l {k)+w 2 {k))) p* cot 5 (s)+p* cot <p d {s) ' 1 ; 

where the quantity F(t) can be read off the Lagrangian in Eq. (34) at tree level 

F{t) = l + D 1 t + 0{t 2 ), t = {k 1 + k 2 ) 2 , 

K = ^(k),| + kV A£=L,(k),|-kY (49) 

Using Eq. (44), we may now perform the integration over the variable Pq in Eq. (48), with the 
result 



F [ k p. T i = p (f) V 6 EnX ° 4p* cot 5{s n ) sm 2 8(s n ) 



(50) 



A Wl (k)w 2 (k)( Wl (k) + w 2 (k) - E n ) (5>(s) + (<P d (s n ))')E n 
On the other hand, the matrix element in Eq. (47) has the following representation: 

(0|0(xo,P,k)£ A (0)|0) = eA(0)^e-^(0|0(0,P,k)|K(P))(K(P)|j(0)|0). (51) 
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Using Eqs. (46), (50) and (51), we get 



\ {En{P)mm = L ^ m g^A _ _ j__ . (52) 

This is the expression of the matrix element in a finite volume. It should be compared with its 
counterpart in the infinite volume, which is obtained by using Watson's theorem: 

(A; 1 ,fc 2 ;out|j(0)|0) = F(t) , \F(t)\ = \F(t) cos 5(s)\ . (53) 

From Eqs. (52) and (53) we finally get: 



|F(t)| 2 = |L 3/2 (K(P)bX0)|0)| 2 ^^ \5'{s) + (0 d ( Sn ))'| 



(54) 



This expression allows one to extract the absolute value of a scalar form factor in the time-like 
region from the measured matrix element (E n (P)\j(0)\0) in a finite volume. Since the phase of 
this form factor, which is determined by Watson's theorem, is also measurable on the lattice, 
we finally conclude that the real and imaginary parts of the form factor can be measured on 
the lattice in the elastic region. 

In the rest frame, the expression in Eq. (54) is similar to the expression obtained in Ref. [21], 
apart from a difference in a kinematical factor which is related to the fact that there a vector 
form factor instead of a scalar one was considered. It can be also shown that, by using our 
method, one exactly reproduces the Luscher-Lellouch formula in moving frames [20,37,38]. 



4 Extraction of resonance matrix elements in a finite 
volume 

Having considered the case of the form factor in the time-like region in great detail, we turn 
to the extraction of the resonance form factor. The part of the Lagrangian that describes the 
interaction with the external scalar field A(x), looks now as follows: 

C A = A(x)j(x) = e 1 A(<S>\$ 1 + ■•■) + e 2 A(<S>\<j> 2 + ••■) + E A(^ 2 ^ 2 + ■•■), (55) 

where ei j2 , E , ■ ■ ■ denote low-energy couplings, and the ellipses stand for the terms with higher 
derivatives. It is seen that, in general, the current j(x) consists of one-body currents and a 
two-body current, whose coupling at lowest order is given by E . 
We make the following choice for the resonance field operators: 

fL/2 

O(x ;P) = / d 3 xe- iPl $ 1 (io ) x)$ 2 (so ) x) ) 

J -L/2 

0%o;Q) = [ L/2 d 3 ye^^l(y ,y)4(y 0l y). (56) 

J -L/2 
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Figure 3: The diagrams describing the quantity M = Mi + M 2 in Eq. (58). Mi and M 2 
correspond to the contributions of the one- and two-body currents, respectively. The wiggly 
line corresponds to the external field A(x). All initial- and final-state interactions are summed 
up in the bubble chains. 

The first-order scattering amplitude of the particles 1 and 2 in the external field A(x) can again 
be calculated using perturbation theory. The pertinent diagrams are depicted in Fig. 3. The 
result of the calculation is (cf. with Eqs. (47,48)): 



(0|(9(x ,P)£ A (0)e%o,Q)|0} = A{0)F(P,Q,x ,yo) , x >0,y <0 



(57) 



where 



F(P,Q,x ,2/o) = ^^ f 

P,q J 



dP dQ c 



e^o M (p 0)P) p. Q 0j Q, q ) e 



-iQoyo 



2ni 2iri 4^!u; 2 (P — i(wi + w 2 ))4:w' l w' 2 (Q — i(w[ + w' 2 )) 



(58) 
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where 



m 2 +[— + pj, w 2 = Jm 2 +\—- 



2 7 z y v 2 

The diagrammatic expansion of the quantity M is shown in Fig. 3. It consists of the contribu- 
tions corresponding to the one-body and two-body currents (see Eq. (55)). Retaining only the 
S-wave contribution in the initial- and final-state rescattering amplitudes, we get: 

M M } • M,. Mi = M[ l) + M[ 2) + M x (3) + M[ 4) , (60) 

where 

M i ] = ( 2w2 ) ri (i + p ' y + q ) i35 f + ( 2w ^ 2 (i - p - y " q ) i3 ^ +P ^ +q ' 
ri|? + P,Q-? + P ) r 2 (^-p,Q-^- P 



M} 2) = — ^ 2 - J S(q*) + — K —, 2 - r J S(q*) + M (2) 

2w 1 (w 1 + w 2 + iQo) 2w 2 {w 2 + wi + iQ ) 

r i p ~i7 + q. ir + q r 2 p -ir-q>ir-q 



M x (3) = V - 2 , 2 J S(p*)+ ± r , , 2 2 x 7 S(p*) + M[ z) 

1 2w' 1 (?I; / 1 +w' 2 + iP ) KF ' 2w' 2 (w' 2 + w[ + iP ) Ky ' 1 

1 ^ 5(p*)(rx(P - 1, Q - 1) + T 2 (P - 1, Q - \))S(q*) 



< } - 73 S 



L 3 ^ 8u;(P - l)w(Q - l)w(l)(u>(Q - 1) + u>(l) + zQ )MP - 1) + + iPo) 

+ M{ 4) . (61) 

In the above expressions, and M x do not contain denominators linear in Qo and Pq, 

respectively, but still include the factors S(q*), S(p*). The quantity M± contains at most one 
energy denominator and both S(q*) and S(p*). These quantities emerge, because initial- and 
final-state rescattering occurs, in general, off the energy shell. Further, 



P \ 2 / / P 



w x = \\m 2 + ( Q- Y + P) ' w 2 = i / m 2 + I Q - — - p 



w[ = x j m 2 + [P-Q + q \ , w 2 = Jm 2 +(-P-Q-q) , (62) 
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p* cot 5(s) + p* cot d (s) ' q* cot <5(s') + <?* cot d '(s') 

(p 2 + p 2 ), ^-(Qg + Q 2 ), 



r 2 * / s ' 2 ^ 27rP ^/ 27r Q 

p = y4 _m ' ^ = \/^ _m ' "r - ' zT' ( ) 

and the I\2 denote the tree-level interaction vertices of the external field A(x) with the fields 

$1,2- 

After projection onto S-waves, the two-body current leads to the following contribution (see 
Fig. 3): 

m P*coU(g) , n . <fcot(s') 



p* cot 5(s) + p* cot d (s) ' ' ' q* cot S(s') + g* cot d '(s') 

+ regular functions in P or Q , , (64) 

where the quantity Z is a low-energy polynomial. 

It can be straightforwardly checked that the sum of all terms in the integrand in Eq. (58) 
do not have singularities at the free two-particle levels. The only singularities are simple poles 
that correspond to the energy levels in the full theory and emerge after the summation of the 
bubble chains. Taking this fact into account and performing the contour integration in the 
variables Pq, Qo by using Cauchy's theorem, we get: 

L T^tl 4 ^i^2^„K + w 2 - E n )(5'{s n ) + (0 d (s n ))') 
327rsin 2 ^( gm ) v ^e^° 

Aw^E^ + w^-E^^s^ + ^s^)')' { ] 

where 

v f P o) J_V iMP-LQ-q + iyp-LQ-i) 

nm[ ' ^ 8w(P - l)w(Q - \)w(\)(w(P - 1) + w(l) - E n )(w(Q - 1) + w(l) - E m ) 

PnCOtS(Sn) rjty „ „ rVl <4 COt S(s m ) . . 

+ —5 — ^^Z(P n ,P;P m ,,Q) — - — — — (66) 



with p* = p*(s = s n ), q* m = q*(s' = s m ). In Eq. (65), the Gottlieb- Rummukainen equation will 
be further used to remove the summations over p, q: 

J_ 1 = Pn cot 5(s n ) 

L 3 ^ 4w 1 w 2 {w 1 + w 2 - E n ) ' 8^^ 

On the other hand, inserting a full set of the eigenstates of the full Hamiltonian, we get 

F(P,Q,x 0?2 /o) = ^(0|O(0;P)|£; n (P))e- B ^(£; n (P)|i(0)|^(Q)) 

n,m 

x e Emy °(E m (Q)\O\0;Q)\0). (68) 
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Further, by using perturbation theory, it is straightforward to show that 



|(0| O (0; P,| £ „,P))| 2 = L> s ™l S $l )y ■ m 

Taking n = m, we readily obtain: 

<^«<»> - £h$Z)y immm K - (P; Q) ■ (70) 

Independently, one may extract the resonance matrix element in the infinite-volume non- 
relativistic EFT by using the procedure described in the introduction. The result is given 
by: 



where 



(P|J(0)|Q> = B R V°°(P; Q) , s, s' s R = 4(m 2 + p 2 R ) , (71) 



8^ B R MlTy/sAPR 



p cot 5 (p) -ip s R -s' 2p R h'(p 2 R ) 



and 



1 1 

h(p 2 ) = p cot 5(p) = — - + - rp 2 H , (73) 

and V°° is obtained from V nm through E n — > P , E m — > Q , s, s' — > and further replacing 
the discrete sum by integration over the variable 1. 

At this stage, we can visualize the problem inherent to the extraction of the resonance ma- 
trix elements. On the lattice, one may measure the quantity (-E , n (P)|j(0)|£' ri (Q)) and extract 
the quantity V nn through Eq. (70). If we were dealing with a stable bound state, in the infi- 
nite volume V nn —> V°°, up to exponentially small corrections. Multiplying with the pertinent 
bound-state renormalization factor, we would directly arrive at the matrix element of the cur- 
rent j(0), sandwiched between the stable bound-state vectors. However, we are dealing with a 
resonance and not with a stable bound state. This means that: 

i) No single E n corresponds to a resonance. We have to formulate a procedure for the analytic 
continuation of the matrix elements into the complex plane. 

ii) The quantity V nn does not have a well-defined limit as L — > oo and E n above the two- 
particle threshold. The 1-loop diagram with an external field, which contributes to the 
M[ 4 \ is the culprit. On the contrary, the contribution from the two-body current, Z, is a 
low-energy polynomial and does not cause any problem. 

In the following sections, we shall explicitly address both of these problems. 
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5 Analytic continuation and fixed points 



In order to avoid kinematical complications, let us first consider the form factor at a zero 
momentum transfer P = Q = 0. The quantity V nn is then a function of a single variable 
p = \J E^/A — m 2 . The questions can be now formulated as follows: 

i) How does one perform the analytic continuation p — > p R in the quantity V nn (p)l 

ii) How does one perform the infinite volume limit L — > oo? 

We shall see below that these two questions are intimately related. 

Let us imagine for a moment that the contribution from the loop diagrams vanishes, so 
that the quantity V nn is given by the two-body current diagram Z only Then, the answers 
to the above equations are trivial. The quantity Z is a polynomial in the variable p 2 : Z = 
Zq + Z\p 2 + • • • . So, one has to first fit the coefficients Zo, Z±, - ■ ■ to the lattice data, and 
then simply substitute p 2 — > p\. The result gives the analytic continuation Z(p 2 ) — > Z{p 2 R ). 
Moreover, since Z(p 2 ) is L- independent, so is Z(p 2 R ), and the final result does not depend on 
the energy level n we started from. 

Let us now see what changes when the one-body current contribution is also included. To 
this end, we first study the analytic continuation of the Luscher equation into the complex 
plane. To ease notation, we restrict ourselves to S-waves and write down the equation (in the 
CM frame) in the following form: 

H P 2) - cot 5(p) = ^Z 00 (1, V 2 ). (74) 

On the real axis, 



p 7T 3 / 2 r] 



V = ^, (75) 

and Eq. (74) determines the energy levels given the scattering phase (or vice versa). Let us now 
look for solutions of this equation for complex values of p. The quantity h(p 2 ) is a low-energy 
polynomial in p 2 , so the analytic continuation is trivial. Furthermore, the function Z 00 (l,ri 2 ) is 
a meromorphic function of the variable rj 2 . Thus, for any given complex value of p, the solutions 
of Eq. (74) determine the trajectories r) n (p), n = 0, 1, ... , in the complex plane (we remind 
the reader that the solutions are not unique). As p — > pr in the p-plane, rj n {p) — > rj n R in the 
//-plane and Eq. (75) becomes a relation that defines L. Our first task is to find all rj n R. 

It is instructive to begin from the 1+1-dimensional case [23]. The counterpart of Eq. (74) 
in this case reads: 

COt S(p) = — COt TTT] . (76) 

The solution of this equation with respect to rj reads: 

% — 1 ~\~ tx 

ri = -— In—— — , x = cot 5(p). (77) 
2ir 1 + ix 
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On the resonance position, we have p — > pr and cot S(p) — > —i. Writing x = — i + e, we get 

V ~ tt" ar S * e - 7T- ln 4" + °( e ) , &sp->PR. (78) 

If we exclude those paths connecting p and pr in the p-plane, which wind around pr infinitely 
many times, then 

lim Rer/(p) < oo , Im rj(p) — > — oo, as p — > pr . (79) 

Recalling the definition of the variable r] (see Eq. (75)), one may interpret the above result (in a 
loose sense) as the equivalence of the mass-shell limit for a resonance (p — > pn) and the infinite- 
volume limit. The same is true for a stable bound state: its energy is volume-independent up to 
exponentially small corrections, so the walls can be safely moved to infinity. Our result shows 
that the same statement holds for a resonance pole (in the 1+1 dimensional case). On the 
contrary, the discrete spectrum above the two-particle threshold is determined by the presence 
of the walls. If one moves the walls to infinity (L — > oo), each given energy level collapses 
toward threshold. The spectrum becomes continuous in this limit. 

What does change in the 3+1 dimensional case? There are so-called finite fixed points with 
\r] n n\ < oo, in addition to the fixed points at infinity which are given by Eq. (79). In order to see 
this, we provide below a numerical solution of Eq. (74) (an analytical solution is not available 
in the 3 + 1-dimensional case). 

The fixed points are the solutions of the equation 

Z o(l;r/ 2 ) + *7r 3 / 2 77 = 0. (80) 
If Im?7 < 0, one may use the following representation of the zeta-function: 

Z o(l;r/ 2 ) = 7r 3 / 2 J-z+ o-Tl e ~^ M } ' n G ^ ' (81) 

^ |n|^0 7r//|n| J 

By using Eq. (81), Eq. (80) can be rewritten as 

6e -2ni V + ^2 + = Q ^ ^ = _L e -2«*»M . ( 8 2) 

In^v^ 

This equation has infinitely many solutions. In order to verify this statement, first assume that 
£(77) = 0. In this approximation, there exists a tower of finite fixed poles parameterized as 

Finally, the equation (82) can be rewritten as 

n + \ ^ln2 % ( y/2 . 



V = 7= 7= In ( 1 + -7— e 2nW2v S(r/) . (84) 

' y/2-1 4tt(V2-1) 2tt(V2-1) V 12 
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Figure 4: Positions of the finite fixed points in the complex ?7-plane. The crosses and circles 
denote 77^ and 7] n n, respectively. The plot is symmetric with respect to the imaginary axis. 



This equation can be easily solved by iteration, starting from 77 = 77^. Note that the series 
for 2(77) contains exponentially suppressed terms and converges very fast. So, truncating the 
sum at some |n| = n m3X can be justified. The numerical solutions indeed exist and are shown 
in Fig. 4. 

Enter the culprit. What are the implications of the above result for the calculation of 
resonance matrix elements? Consider a simplified expression for V nn (0, 0) in Eq. (66), setting 
I\ = 1, r 2 = and Z = (the low-energy polynomials in the numerator do not alter the 
analytic properties of the diagram we are interested in, and the part containing Z is trivial and 
was considered already). All we have to consider is the expression 



Kn(0,0) 



-y 

L 3 ^ 



(2t/7(1) ) 3 (2w(l) 



m 



p 



8tt£ 3 p 2 



pcot 0(77) 



2 2 

m — p 
8nE 3 p 2 



pcot S(p) 



E n ) 2 



1 

32nEp 
1 



d 

Je 



p Z m {\\rf 



8nE 2 
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3/2 



E—E„ 



ZlixEp 



;i + cot 2 0(r/))7/0 / (7/) 
+ cot 2 5(p)) rjfp'irf) 



+ ■ ■ 



E—E n 



+ 



(85) 



E — En 



where 4>(rj) = (j) d=0 (E 2 ), E = 2^J m? + p 2 , and the ellipses stand for the terms which vanish 
exponentially with L. In the last line of Eq. (85), Luscher's equation cot 0(77) = — cot 5(p) 
was used. The quantity V nn (0, 0) is a function of the variable p, so one can write V^O, 0) = 
V nn {p). It is now legitimate to ask how the analytic continuation of the above expression in p 
is performed and what is the result of this continuation. The expression in Eq. (85) consists of 
two terms. It can be verified directly that the first term is a low-energy polynomial in p 2 (up 
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Figure 5: Two different paths in the //-plane and corresponding paths in the p-plane. It is seen 
that the paths converge to the same point pr in the p-plane. 

to a trivial overall factor p~ 2 ). The analytic continuation of this term is straightforward and 
leads to 

2 2 2 2 

m — p . . m — p R 

———pcot5(p)^-i—^, ^P^Pr. (86) 
m & p 8ns R Pr 

It is easy to check that this result exactly coincides with the result for the loop diagram 
calculated in the infinite volume (i.e., replacing summation by integration in Eq. (85)), on the 
second sheet. Consequently, if the second term, continued to p — pr, vanishes, the analytic 
continuation of the whole vertex diagram to the pole on the second sheet will yield the same 
vertex evaluated in the infinite volume. This would be the statement that we are after. 

Let us assume for a moment that it is possible to find a procedure to perform such an 
analytic continuation in the second term of Eq. (85). We choose some path in the complex 
p-plane approaching the pole at p — p R . Suppose first that, moving along this path, the 
variable r\ = r](p) approaches the infinite fixed point |Re^ n ij| < oo, lmrj n R — > — oo. Using the 
representation for the zeta- function given in Eq. (81), it can be easily checked that the second 
term in Eq. (85) indeed vanishes if 7/ tends to the infinite fixed point. 
Imagine now a path that ends at a finite fixed point. Parameterizing this path as 

cot <f>{ri) +i = \( v - Vur) + 0((t] - r] nR ) 2 ) , (87) 

where A is a finite complex constant. Then, in the vicinity of the fixed point, 

<j>{rj) oc ln(r/ - rj nR ) + const + 0(r) - r] nR ) , (j)'{rj) oc h const + 0(r) - r) nR ) . (88) 

V — VnR 

From the above equations it is evident that the product 7/(1 + cot 2 (f)(r)))(f)'(r)), rather than 
vanishing, tends to a constant at the finite fixed point. In other words, if during the analytic 
continuation, the variable 7/ gets caught by a finite fixed point, the result of the analytic 
continuation is different from the vertex function in the infinite volume and one is in trouble. 

In order to understand this result better, let us consider some point 7/ on the real axis and 
two paths, connecting r] to an infinite and to a finite fixed points, respectively (see Fig. 5). 
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These paths are mirrored by pertinent paths in the p-plane. Since we have assumed that 
there is only one resonance pole at p = pr, both paths in the p-plane start at the same point 
p = Po corresponding to rj = r/o and end at the same point p = Pr. The result of the analytic 
continuation is, however, different along these paths, rendering an unambiguous determination 
of the vertex function at p = pr impossible. 

The problem, which was discussed above, looks complicated but has a particularly simple 
solution. Let us go back to the last line of Eq. (85). It is immediately seen that the ambiguity 
is caused by the expression r)(f)'(r)), which is contained in the second term and depends on the 
energy level index n. Moreover, the form of this expression is universal (it does not depend on 
the interaction). Consequently, measuring the vertex function for two different energy levels n 
and m, and forming the linear combination, 

y,x = Vnn(p)Vm(j)'(r)m) ~ Vmm(p)Vn<j>' (Vn) , gg s 
Vm4>'(Vm) ~ Vn&'iVn) 

one may immediately ensure that the culprit disappears. Namely, V(p) is a low-energy poly- 
nomial in the variable p 2 up to a factor p~ 2 , it does not depend on the energy level (up to 
exponentially suppressed contributions), and its analytic continuation p — > pr into the complex 
p-plane yields the infinite-volume vertex function. To conclude, the problem with the analytic 
continuation was circumvented by measuring the matrix elements for two different energy levels. 

Finally, we would like to note that the problem is milder in the case of 1+1 dimensions, 
see Ref. [23]. First, there are no finite fixed points and no ambiguity emerges. Second, in 
Ref. [23] it has been shown that the problematic contributions in 1+1 dimensions can be fit 
by a polynomial in p (not p 2 ) with n-dependent coefficients, so the analytic continuation still 
can be performed (although it is a more subtle affair now, see Ref. [23] for the details). No 
similar statement exists in the case of 3+1 dimensions. The subtraction trick can be used in 
1 + 1 dimensions as well, making the fit more straightforward (at the cost of measuring two 
energy levels instead of one). 



6 Matrix elements at nonzero momentum transfer 

We finally turn to the resonance matrix elements for non-zero momentum transfer. It is con- 
venient to work in the Breit frame P = — Q. The vertex function in the infinite volume using 
dimensional regularization is given by (we again neglect the numerators which do not affect the 
analytic properties): 




(2ir) d 8w(L>(P - l)w(P + 1) (w(P - 1) + w(l) - P )(w(P + 1) + w(l) - P ) ' 



p* = y^-™ 2 , s = P 2 -P 2 , (90) 

where the Po —> Po+iO prescription is implicit. The finite- volume counterpart of this expression 
contains a sum over the discrete momenta 1 instead of an integral. We note here once more 
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that a particular prescription is used to calculate this integral: the integrand is first expanded 
in powers of the momenta, integrated over 1 and the resulting series is summed up again. 
Using this prescription, one may present the above integral in the following form (consult, e.g., 
Ref. [32] for the technical details of similar calculations): 

/ = h + h, 

1 f d d \ 1 

h 



2P J (2*)<4PH ( ,_ ?)a 1 (P(1 _ f))! _ (p>)2 



p2 



1 



a + f) 2 -4( p ( 1 + !)) 2 -(^) 2 

d d \ f 1 



2P J (2n) d \ 2w(P + l)(w(P - 1) + w(P + 1)) (] _ E)2 _ _L (p(1 _ £ ))2 _ (r)2 
1 1 



p2 



2 W (P - l)(w(P - 1) + w(P + 1)) (1 + |)2 _ 1 (p(1 + |))2 _ (p , )2 



(91) 



p2 
M) 



Explicit calculations yield the following result (on the second sheet): 

1 , -«'Vi|P| 

arctan ■ 



16ttP |P| 2PoP* 



32iry/s 7_! * \ io 2 (wi + «7 2 ) w'l {w'{ + O 



where 



= + y )' f i + g^u 2P °'y » (i ± n + (i ± n v> 1/2 . m 
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Now let us consider the same quantities in a finite volume: 
2*>»^ m \ _ 1 _ 



p2 



1 



1 



(l+f) 2 -^2(P(l+f)) 2 -b*) 2 



7 l _ _J_J_W 

J 2 - oprsL 



2P L3 ^ 1 2^(P + 1)(«;(P - 1) + u,(P + 1)) (1 _ P)2 _ i_ (P(1 _ p ))2 _ (r)2 

Pn 



(94) 



2»(P - 1)(„(P - 1) + u>(P + 1)) ( , + ?)a _ 1 (p( , + p ))2 _ (p>) 



Neglecting partial-wave mixing in the finite volume, the quantity l£ can be rewritten as 

327rys y_i I ^ 2 (^i + w 2) w'{(w 1 + w 2 ) J 

Using the Gottlieb-Rummukainen equation, it is straightforward to ensure that 1% is a low- 
energy polynomial and its analytic continuation to p = p^ gives the infinite- volume result Ii- 
On the contrary, I± does not have the same property. For this term, we use the following trick. 
We define: 

It = if + (If - If, = if + s , If = (M arctan oo t , w . (96) 

The quantity If is a low-energy polynomial (up to a trivial overall factor (p*)~ 2 ), and its 
analytic continuation to the pole on the second sheet gives I±, which is the value of the integral 
in the infinite volume. Further, the quantity g is dependent on the energy level, and is universal 
(all derivative interactions factor out). Consequently, measuring the vertex function for two 
different energy levels n and m in the Breit frame, and forming the linear combination 

y,*s = Vnn(p*)g m (p*) ~ V mm (P*)9n(p*) ^ 

g m (p*) - g n (p*) 

one sees that the culprit cancels out: V(p*) is a polynomial up to a factor (p*)~ 2 , and its 
analytic continuation to the resonance pole yields the vertex function in the infinite volume. 



7 Conclusions 

i) In this paper, by using the technique of the non-relativist ic effective Lagrangians in a 
finite volume, we were able to formulate a procedure for extracting the resonance matrix 
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elements on the lattice. The derivation was restricted to the case of isolated resonances, 
lying in the region of the applicability of the effective-range expansion. 

) As a demonstration of the usefulness of the non-relativistic EFT approach, we have re- 
derived the Liischer equation in the moving frame (Gottlieb- Rummukainen equation), as 
well as the relation of the time-like form factor to the matrix elements measured on a 
Euclidean lattice. 

) A resonance pole is extracted in the following manner: by performing the measurement 
of the energy levels at different volumes, and using Liischer's formula, one extracts the 
function p cot S(s) at different values of p. In the region of applicability of the effective- 
range expansion, which we have assumed here, this function is a polynomial in the variable 
p 2 : pcot5(s) = Aq + Aip 2 + • • • (for simplicity, we consider the S-wave). The fit to the 
lattice data determines the coefficients Aq,Ax,---. The resonance pole position is then 
determined from the equation 

p R cot 5 (s R ) = A + Axp 2 R H = -ip R . (98) 

Note that a shortcut version of this procedure is to determine the zero of the function 
pcotd(s) and to relate the width of a resonance to the derivative of this function. At 
present, this shortcut version is routinely used to study the resonance properties on the 
lattice. For narrow resonances, both procedures give the same result. 

) The case of the resonance form factors is more subtle. It has been demonstrated that a 
straightforward analytic continuation of the matrix elements of the current between the 
eigenstates of the Hamiltonian in a finite volume does not allow one to determine resonance 
matrix elements unambiguously in 3+1 dimensions, and the infinite volume limit can not 
be performed. 

) The way to circumvent the above problem is to measure the matrix elements for two (or, 
eventually, more) eigenstates. The extraction of the matrix element proceeds in several 
steps: 

— Use the Breit frame, then extract matrix elements between at least two different 
eigenstates, labeled by n, m, by using Eq. (8) (or its counterpart for excited states). 

— Using Eq. (70), extract the quantities V nn (p*),V mm (p*) with p* = {{E 2 - P 2 )/4 - 

1 /2 

m 2 ) , and E = E n or E m . Note that, in the Breit frame, V nn (p*),V mm (p*) depend 
only on p*, as P is fixed. 

— Form the linear combination V(p*), using Eq. (97). Fit the results of the measurements 
for different values of L by using the formula 

W = + D * + + • • ■ • (") 

— Calculate V°° = V(pr) by simply substituting p* = pr in the above expression. 

— Finally, calculate the resonance form factor in the infinite volume by using Eq. (71). 
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vi) The procedure described above demands that the matrix elements between the eigen- 
states are measured on the lattice at several different volumes and at least for two 
different eigenstates. We realize that, at present, this requirement is rather chal- 
lenging. However, in our opinion, it is still important to have a clearly defined and 
mathematically rigorous procedure, which will allow for a clean extraction of resonance 
form factors in the future. Turning the argument around, our discussions demonstrate 
that the existing lattice results for the resonance matrix elements should be put under 
renewed scrutiny. 

vii) It would be interesting to extend the discussion to the case of twisted boundary 
conditions, which have proved advantageous in the calculations of form factors. Non- 
relativistic EFT is ideally suited for this purpose. We plan to investigate this issue in 
the future. 

viii) In this paper, one has assumed that the effective-range expansion is valid for the 
energies where the resonance is located. It would be interesting to extend the range 
of applicability of the approach, by using e.g. conformal mapping. 
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